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ABSTRACT 

General relativistic magnetohydrodynamic (GRMHD) flows along magnetic fields threading a black 
hole can be divided into inflow and outflow parts, according to the result of the competition between 
the black hole gravity and magneto-centrifugal forces along the field line. Here we present the first 
self-consistent, semi-analytical solution for a cold, Poynting flux-dominated (PFD) GRMHD flow, 
which passes all four critical (inner and outer, Alfven and fast magnetosonic) points along a parabolic 
streamline. By assuming that the dominating (electromagnetic) component of the energy flux per flux 
tube is conserved at the surface where the inflow and outflow are separated, the outflow part of the 
solution can be constrained by the inflow part. The semi-analytical method can provide fiducial and 
complementary solutions for GRMHD simulations around the rotating black hole, given that the black 
hole spin, global streamline, and magnetizaion (i.e., a mass loading at the inflow/outflow separation) 
are prescribed. For reference, we demonstrate a self-consistent result with the work by McKinney in 
a quantitative level. 

Subject headings: galaxies: active-galaxies: jets — Magnetohydrodynamics (MHD) — Black hole 

physics 


1. INTRODUCTION 


Relativistic jets emerging from accreting black hole 
systems have been observed in active galactic nuclei 
(AGN), micro-quasars (stellar mass black hole X-ray 
binaries), and presumably gamma-ray bursts (GRBs). 


AGNs are > io — 20 (IJorstad et al. 

20051 Cohen et al.l 

2007; Gu 2009; iPushkarev et a,1112009 

Lister et al.l 2013), 


and the values could be h igher for some blaz ars (see 
iGonal-Krishna et all 120061 iHovatta et all 120091) . The 
Lorentz factors of mic ro-quasar jets are l ower, mostly 
with T ~ 2 — 10 (e.g. iFender et all 120041 : 1 Corbel et al.l 
20021). but still there are a few found to have T > 10 (see 
Miller-Jones et aI1l2006lf . Jets in gamma-ray bursters are 
supposed to be ultra-relativistic, and their Lorentz fac¬ 
tors can be as high a s ~ 100 — 1000 (see, e.g. iLvubarskvl 
120101: iLvutikovl 120111 ). How jets become relativistic af¬ 
ter being launched from nearby black holes is a long¬ 
standing issue. Electromagnetic or magnetohydrody¬ 
namic (MHD) mechanisms are frequently invoked to ex¬ 
tract energy and mo ment um from the black hole and 
accretion disk (e.g.. lMeier et al.ll200lL for reviews). One 
of the key issues to be addressed is the potential of the 
MHD flow acceleration up to a high bulk Lorentz factor 
of F > 10. 


An ideal engine to power relativistic jets is a spinning 
black hole. Close to the black hole, the rapid winding 
of the azimuthal component in large-scale magnetic 
fields due to the frame-dragging inside the black 
hole ergosphere results in a counter torque (induced 
by the Lorentz force) against a black hole rotation. 
The energy that the black hole spent to perturb the 


field line can be propagated outward in the form 
of torsional Alfven waves, thus extracting the b lack 
hole energy electromagnetically dBlandford fe Znaiekl 
H977l hereafter BZ77). However, because the en¬ 
vironment around an accreting black hole is not a 
perfect vacuum (contrary to the force-free treatment 
in BZ77), the general relativistic magnetohydrody¬ 
namics (GRMHD), which consist of electromagnetic 
and fluid components, provide a more general pic- 


jets in both theoretical approaches (e.g. Camenzind 

1980a] 

U 1987]; iTakahashi et all 

L99Q; [Fendt & Greiner 

200j 

Fendt & Camenzindl 1996; 

Fendt & Ouyed 200-1) 


2000 

: IMizuno et al.l 2004 IMcKinnev & Gammie 

2004 

: lHawlev & Krolik 20061: McKinnev 

2006; 

Beckwith et al. 20081 Tcheklrovskov et al. 20101, 

2011; 

Tcheklrovskov & McKinnev! 2012]). 


The overall configuration of an accreting black hole 
system is schematically illustrated in Figure |T| (see 
also G RMHD simulatio n s for magnetized acc r etion , 
e.g. iMcKin nevI 1 2001)1 : IMcKinnev fe; Gammiel 12004 
lHawlev fc KroliU 1200(1) . Ordered, parabolic lines are 
developed near the funnel, which is confined by the 
corona and/or accretion flow. Due to the relative ab¬ 
sence of accreting materials, the funnel region is Poynt¬ 
ing flux-dominated (PFD). The fluid loading onto the 
field is accelerated inward (or outward) if the black hole 
gravity force is larger ( or smaller) than the “m agneto- 
centrifugal” forces (e.g. iSadowski fc Sikoral 120101 in the 
case of the accretion disk). As poin ted out in the theoret¬ 
ical work of iTakahashi et al. ( 1990 ). black hole rotational 
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energy can be extracted outward by a PFD GRMHD 
inflow. This is a direct result of that the electromag¬ 
netic components dominating the GRMHD flow, and 
the electromagnetic component is responsible for extract¬ 
ing the black hole energy, similar to the BZ77 process. 
The outward energy flux, after being extracted from the 
black hole, is expected to propagate continuously out¬ 
ward throughout the magnetosphere from the inflow re¬ 
gion to the outflow one. In this paper, we focus on the 
PFD GRMHD flow in the funnel region, including both 
the inflow and outflow parts. 

For comparison, let us quickly consider the case when 
the GRMHD flow becomes fluid-dominated. In that 
case, the energy flux is dominated by the fluid compo¬ 
nent, and therefore it has an inward direction for inflow, 
but outward for outflow (c.f., the energy flux direction 
shown in FigureQ]). Such discontinuity of the energy and 
momentum fluxes implies that the outflow is accretion- 
powered, which is constrained by the energy input from 
the disk/corona. The switch-on and switch-off of the 
extraction of the black hole energy (inflow) may closely 
relate to the launching and quenching of rela tivist ic jets 
(outflow) (e.g. lPu et al.ll2012t iGlobus fc Levinsonll201.Hl. 

Prio r to the GMRHD studies mentioned. IPhinnevI 
(11983H considered the inflow and outflow along a 
monopole field jointly by the conservation of the to¬ 
tal energy flux per flux tube. In this pioneering work, 
they consider energy extraction from the black hole via 
BZ77 process (the inflow part), and the Michel’s “mini¬ 
mum torque solution” (lMichelll969 ). in which the fast(- 
magnetosonic) point is located at infinity (the outflow 
part). We, however, suggest that a more realistic situ¬ 
ation can be considered: the black hole energy extrac¬ 
tion process in the framework of GRMHD, and a type of 
parabolic GRMHD flows as a result of external pressure 
confinements provided by the corona/accretion. Recent 
observational evidence also supports this idea; nearby ac¬ 
tive radio galaxy, M87, exhibits the parab olic, stre amline 
up t o ~ 10 5 Schwarzschild radius (lAsada fe Nakamural 

i2oTl . 

Furthermore, we are interested in the case that the 
fast point of the outflow is located at a finite distance. 
This consideration is directly related the conversion from 
Poynting to kinetic energy fluxes of the flow and there¬ 
fore the jet acceleration. Poloidal magnetic flux is re¬ 
quired to diverge sufficiently rapidly in order for most 
of the Poynting flux to be converted into the kinetic en¬ 
ergy flux beyond the fa st point (also known as the mag¬ 
netic nozzle effect (e . g. iCamenzin dl [19891: ILi et al.lll99i 
iBegelman fe Lilll994 iTakahashi fc Shibatall 19981) . 

iBeskin fe Nokhrinal (l2006h examine the acceleration of 
the jet along a parabolic streamline by introducing a 
small perturbation into the force-free field. As a result, 
the fast point is located at a finite distance. This indi¬ 
cates how plasma loading in the flow plays a role in ac¬ 
celerating the flow, as well as a conversion from Poynting 
to kinetic/particle energies. They consider the behavior 
of the outflow in the flat spacetime. However, we are 
interested in both the inflow and outflow near a black 
hole. 

All of these theoretical works provide important pieces 
toward a picture that includes the following process along 
the field line: (i) in the inflow region the rotational energy 
of the black hole is extracted outward by the GRMHD in¬ 



Fig. 1. — Schematic illustration of a Poynting flux-dominated 
(PFD) GRMHD flow confined by the accretion flow and its corona. 
The outward-streaming curves indicates ordered, large-scale mag¬ 
netic fields that thread the black hole event horizon. The inflows 
and the outflows (represented by thick white arrows) are along 
the field lines, and they are separated by the separation surface 
(marked by a dashed line). The energy flux (represented by a gray 
arrow) is outward in both the inflow and outflow regions, as the 
black hole rotational energy is extracted and transported outward. 
The static limit (dashed curve) and the light surface (solid curve) 
outside the black hole (black region) are also shown. 


flow, (ii) at the the inflow/outflow separation surface the 
extracted energy flux is carried out continuously, and (iii) 
in the outflow region the flow passes the fast point, and 
hence the bulk Lorentz factor increases. Although this 


picture has been already recognised in the quasi-steady 
state in GRMHD simulations ie.g. livl cKinney] | 200fil 
iMcKinnev fc Gammiel 120041 : lHawlev fc Krolikl 120061 ). no 


steady solution is available in the literature. 

In this paper, we present the first semi-analytical work. 
We consider the energy extraction from the black hole 
via the GRMHD (inflow) , and t h e pertur bed force-free 
parabolic field line in IBeskin fc Nokhrinal (1200611 (out¬ 
flow). With given black hole spin, field angular velocity, 
and magnetization at the separation surface, we are able 
to to constrain the outflow solution by the inflow solu¬ 
tion. For a reference, we adopt simila r parameters re- 
ported in the GRMHD simulation of iMcKinnevI (120061 
;hereafter M06). Our semi-analytical solution passes 
all the critical points (inner and outer, Alfven and fast 
points), and agrees with the inflow and outflow proper¬ 
ties along a mid-level field line in M06. 

The paper is organized as follows. In §2, we outline 
the GRMHD formulation and the wind equation (WE). 
In §3, with the consideration of the conservation of en- 
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ergy flux in inflow and outflow region near the separation 
surface, we discuss the matching condition to connect the 
inflow and outflow part of a PFD GRMHD flow. In §4, 
we introduce our model setup. We adopt similar param¬ 
eters to those reported by M06 and compare the solu¬ 
tion obtained by the matching condition with that of the 
time-averaged GRMHD numerical simulation results in 
M06. Finally, a summary is given in §5. 


2. STATIONARY, AXISYMMETRIC MHD FLOW IN KERR 
SPACETIME 

2.1. Basic Formulae 

The theory about stationary and axisymmet- 
ric ideal G R MHD f lows has be en in several works 
(iCamenzin dl Il986allbl . 119871: iFendt fe Camenzindl 
19961: Fendt fc Ouved 12004 iTakahashi et al.l 1 1 9901 
Fendt fc Greiner! 12001 1. For completeness, in this 
section we summary and present the necessary formulae 
for the purpose of this paper. 

The natural unit system is used throughout this work. 
As c = G = M = 1, the gravitational radius r g = 
GM/c 2 = 1, where c is the speed of light, G is the grav¬ 
itational constant, and M is the mass of the black hole 
(conversions from the c.g.s. units to the natural units for 
the ph ysica l variables here can be found in tables 3 and 
4 in lPu et all (120121 11. The flows occur in a background 
Kerr space-time, which is stationary and axisymmetric. 
For a metric signature [— + + +], the Kerr metric (in 
Boyer-Lindquist coordinates) reads 

A — a 2 sin 2 9 , 9 4ar sin 2 0 
ds 2 = --- dt 2 --- dtd(f> 


+ Asmd_ ^ + S dr2 + E dQ 2 # (1) 

2.J Za 

where a = J is the angular momentum of the black hole, 
A = r 2 — 2r + a 2 , E = r 2 + a 2 cos 2 9, and A = (r 2 — 
a 2 ) 2 + A a 2 sin 2 9. 

We also assume that the flow is cold. For a highly- 
relativistic flow, the thermal pressure p is insignificant 
compared with the rest-mass energy density and the ki¬ 
netic energy density in the fluid, and hence the cold limit 
is justified. 

The flow is magnetized and the stress-energy tensor of 
the fluid has two components: 

T^ = T^ + T^, (2) 

where the fluid component is given by 

= puV , (3) 

and the electromagnetic component by 

t Za = ^ - \g^F°PF a ^ , ( 4 ) 


where u^ is the 4-velocity of the fluid and p is the rest- 
mass energy density. The electromagnetic field tensor 
Fy V satisfies Maxwell’s equations, and the proper num¬ 
ber density n (= p/m, where m is the rest-mass of the 
particles) satisfies the mass continuity equation. 

Under the ideal MHD condition, a stationary and ax¬ 
isymmetric flow obeys four conservation laws: 

(mA) ;M = 0, (5) 


£ 

II 

o 

(6) 

II 

o 

(7) 

(g^). v = o , 

(8) 

where = d/dt and = d/d(f> are 

the Killing vec- 


tors. These conservation equations (equations [5] 0) give 
four conserved quantities along a streamline. By denot¬ 
ing the poloidal stream function as \H(r, 0), they are: (i) 
the angular velocity of the field line, Dp (’P), (ii) the par¬ 
ticle number flux per unit magnetic flux (mass loading), 
77(H/), (iii) the total energy of the flow per particle, E{\ k), 
and (iv) the total angular momentum per particle , L('F) 
(iCamenzindll 1 986al lbf. I 1 9871 ITakahashi et alJfl990t) : 




Ftr 

F r p 


F t e 
F h ’ 


»7W 


y/—g nu r V~~9 nu6 

Fop F r p 

\J—g nu 4 (D — Sip) 

FVe ’ 


(9) 


( 10 ) 


E(^)=Efl + E E m 

=-pu t - -— V~g f 

47T?7 


r0 


= -gu t 






( 11 ) 


L('F) = Lfl + Lem 


Bp 

= W "4^’ 

(12) 

with y/—g = Esin A Here, 12 = u* /it 4 is the fluid angular 
velocity and p is the relativistic specific enthalpy, which 
becomes m in the cold limit. The covariant magnetic 

field observed by a distant observer with u v = 
is given by 

By, = ^euy a pF ap u 1 ' , 

(1,0,0,0) 

(13) 

and its toroidal component is given by 


Bp = V=g F re = (A/E) sin 9F rB , 

(14) 


where e U y a p = \J~-g [vpaj3\ is the Levi-Civita tensor, 
and [i/pa/3\ is the completely antisymmetric symbol (see 
Appendix). 

The outward energy flux in the flow is 

£ r = —T[ = nEu r , (15) 

and the outward angular momentum flux is 

C r = -1% = nLu r . (16) 

Splitting them into fluid (i.e. Eee and Lee) and the elec¬ 
tromagnetic components (i.e. Lem and Lem) gives 

£ r = £ F L+ W 

= nKpLM r ’ + nEEMU r 
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and 


= —nputu r 


B$ 

4n Esin# 6 


£ r = £ FL r + £ EI / 


= niFLM r + nLEMU r 


= —npu ( j ) u r 


cr 

6 EM 

ftp 


(17) 


(18) 


(iTakahaslii et al.lll990f ) . where 


Kq — f + ^gtip^F + 9tt) 

T x 2 

K 2 = ( l-fi F - 


E 


K 4 = 


-1 


9t(j> 9tt9<t><t> 


0 L 

9<t><l> + 2 g t< j,— + gtt 



(25) 

(26) 

.(27) 


As initially proposed bv lTakahashi et all (Il990f ). the case 
of £ r > 0 for inflow (which requires a negative total en¬ 
ergy) is known as the MHD Penrose process. For later 
studies, the term is instead used to indicate a nega- 
tive energy orbit of t he fluid componen t , ££ T > 0 (e.g. 
Hirotani et al.l 1 19921 IKoide et all 120021 IS ernenov et ali 


2004 iKomissarovI 20051: IKoide fe Baball2014l)~ 

The bulk Lorentz factor of the flow for a distinct ob¬ 
server can be defined by 


r = V-9tt 


(19) 


If all the energy in the Ponyting flux is converted to 
the fluid’s bulk (kinetic) energy at a large distance, the 
terminal Lorentz factor will be 


r 


oo 


E 

9 


£ r 

pu r 


( 20 ) 


In addition, the angular velocity of the fluid at a large 
distance will be 


u 


<$> 

OO 


L 

9 


a_ 

pu r 


( 21 ) 


Equations © and (1211) therefore provide the upper limit 
of the terminal Lorentz factor and angular velocity of the 
fluid at large distances. 


2.2. Wind Equation 

The streamline of the flow is represented by the func¬ 
tion \l/(r, 6) = const. The WE (i.e., the relativistic 
Bernoulli equation), describing the fluid motion along 
the streamlines can be obtained using the normalization 
condition u a u a = — 1. The WE therefore has the form: 

U 2 + 1=(J) U g (r,0), (22) 

where the poloidal component of the 4-velocity is given 

by 

Up = u-’uj , (23) 


The Alfven Mach number M\ is given by 
Ml = 4t rprf/n = 4npn = 4npp (^-^ , (28) 

where the re-scaled poloidal field is 


B p = B J Bj ( g tt + gt</,Q f) 
1 


A sin 9 


[g^r + s'We] 


(29) 


Along the streamline several characteristic surfaces can 
be defined. Their definition and properties are summa¬ 
rized in the Appendix. 

The conserved quantities E, L and p can be expressed 
in terms of three system parameters: (i) the launching 
point of the flow, r*, (ii) the location of the Alfven sur¬ 
face, fa, and (iii) the magnetization parameter at the 
launching point, cr*. Explicitly, the relations are 


L 

E 

E 2 


9t<t> + gu^F 

gtt + gtp^F 


r A 



<7 


★ 




2 

★ 


47rm (nup^/—g) | Tt 


( I>* 1 

47tto \p\ ’ 


(30) 

(31) 

(32) 


where E = E/p, the flux function $ = B p ^/—g, and $* 
denotes that it is evaluated at r = r*. In terms of these 
parameters, the Mach number can be written as 


M a 


-Up 


cr* 


| 1/2 , 


(33) 


where / = 4>*/4>. Note that u p = 0 at r = r * has been 
assumed in deriving the relation (1311) . In addition, the 
relation (1321) implies that knowing the mass loading p is 
equivalent to knowing the a*. 

In the cold limit, the WE is a polynomial equation of 
4-th order in u p : 


with the summation over the poloidal indices j = {r, 6}. 
The term U g (r, 9) in the right-hand side of equation 
(l22l) . which is evaluated along the magnetic field line 
in the calculation, is related to the conserved quanti¬ 
ties, and its explicit express ion depends on the as sume d 
background space-time (seelCamenzindl (fl986albl. Il987l ): 
IFendt fe Camenzindl (11996H : lFendt fe Ouvedl ( 2004T) for 
the Minkowski and S chwarzschild spacetime s; and 
iTakahashi et al.l (I1990H : IFendt fe Greineil (|200lH for the 
Kerr space-time). 

In a Kerr space-time, we obtain 


U g (r, 9) 


K 0 K 2 - 2 K 2 M\ - KiMi 
{Ml - K 0 ) 2 


(24) 


^A«p = 0. (34) 

i=0 

The coefficients Ai(r; T, $, r*,r\, cr*, Hf) are given by 

Ai = C 2 , 

A3 = —2C 2 Rq , 

a 2 =cI + k 2 + e 2 k 4 cI 

Al = -2C 2 K 0 + 2E 2 K 2 C 2 

a 0 = k 2 -e 2 k 0 k 2 , 

where C 2 = yj^^g f /cr* (IFendt fe Greineill20011) . 
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3. MATCHING CONDITION OF THE INFLOW AND 
OUTFLOW 


In the work of IPhinnevl (119831) . the matching of the 
inflow and outflow parts of the flow is constrained by the 
conservation of the energy flux per magnetic flux in the 
inflow and outflow region 


(^7-E) inflow — (l?^)outflow 


(35) 


Remember that both ry and E of the inflow and outflow 
are constant. 

Consider equation (1351) at the separation surface, r s , 
for PFD flow (E Eem O Efl), we further consider 


(^)inflow (Rem) — (^outflow (l^EM ) 


(36) 


to be the matching condition of the inflow and outflow. 
The superscripts ” (or “+”), respectively, denote the 
physical value computed at the location very close to r s 
in the inflow (or outflow) region, that is, r —> r~ (or 
r —> r+). After some algebra, equation (l36l) can also be 
expressed as 

(^em) = (£em) + i (37) 


or 




(38) 


Equation (1371) implies that the matching condition we 
adopt is equivalent to the statement: the outward Poyt- 
ing energy flux is continuous at the separation surfaccQ. 
Equation (l38l) revels that such condition guarantees that 
the toroidal field is continuous at the separation point, 
provided that I2 f is the same constant in the inflow and 
outflow region. 

It is interesting to note that the matching condition 
does not require that ry or Eem should be continuous 
when crossing r = r s . That is, if we define 


6 = 


(^/) inflow 


(Eem) + 


(ct*) + 

(^/)outflow 


(Eem) - 


(a*) - 


(39) 


6 is not necessary for unity. The last relation in equation 
(1591) is obtained with the help of equation (15121) . Never¬ 
theless, due to following reason, the outflow can still be 
properly constrained by the inflow even with the uncer¬ 
tainty of 5. 

Consider a flow along a prescribed, hole-threading 
poloidal_field line with some specific angular velocity 
field. iZnaiekl dl977t) showed that, due to the regular¬ 
ity requirement at the event horizon, r + = 1 + Vl — a 2 , 
the derivative of the stream function T is finite and B$ 
satisfies 


B cf ,(r + ,6) = — sin$- 


' + 


a 2 cos 2 6 


(fl H — , (40) 


where f 1h is the angular velocity of the hole. As a 
result, (E<y,) _ is insensitive to different value of (cr*)“ 
((E^) - ~ const.). From dynamical point of view, this 
can be understood by the fact that the fast point of a 
PFD GRMHD inflow is always located close to the black 
hole event horizon (Appendix). 

For outflow, however, there is no constraint at infinity, 
and therefore {B^) + depends on (ct*) + more strongly. 

1 cf. equation (1351) gives (£ r )~ = (£ r ) + • 


Again, from dynamical point of view, the relatively 
strong dependence can be understood by the fact that 
the fast point of the outflow can vary from finite dis¬ 
tance to infinity. Because the uncertainty of the S is 
introduced by the uncertainty of (cr*) - , instead of (cr*) + 
(see also §4.2), the outflow can still be well constrained. 
The matching condition then play the role to constrains 
the outflow by singling out the outflow solution that sat¬ 
isfies (E^,) - = (E^,) + . 

4. FLOW ALONG A PARABOLIC FIELD LINE WITH A 
FINITE-DISTANCE FAST POINT 

4.1. Model Setup 

In general, the field configuration should be con¬ 
sistently determined by solving the trans-field equa¬ 
tion, (i.e. the Grad-Shafranov equation). The trans¬ 
field equation in cold limit involves the stream func¬ 
tion T, and the derivative of the amsrawesquanti- 
ties. dti.E/d'fy, d/n/d^, dE/dip, dL/dyb (|Nitta et alJll991b 
iBeskin &; Par’evlll993l) . However, solving the trans-field 
equation analytically is very challenging and beyond the 
scope of this paper. 

On the other hand, we are interested in the case where 
the fast point of the outflow is located at a finite dis¬ 
tance. It is therefore essential to consider an additional 
modification on the original force-free field line due to 
the MHD flow. We leave a better consideration of field 
configuration f or a future work, and ad opt the streamline 
function in IBeskin fc Nokhrina (12006? ) as the prescribed 
parabolic field 

T = T 0 + ef , (41) 

where To is the the flat spacetime parabolic force-free 
field generated by the toroidal surface current distri¬ 
bution, I = C/47rr(l + Hpr 2 ) 1 / 2 , on equatorial plane 
(iBlandfordll 19761 iLee fe Parkll2004f) . 

nC 

To = sinh - (fip^l — cost))) , (42) 


and 


ef = e TTCflpr sin 6 , e<l 


(43) 


is the perturbation introduced by the MHD effect. The 
constant C is assumed to be unity. Note that by the 
help of the relation sinh -1 (:r) = ln(ir + V x 2 + 1), To 
is proportional to r( 1 — cos6 ), which is the same as 
the dominating term of the parabolic fielcQ in BZ77. 
In addition, V • B = 0 is guaranteed. It is shown in 
IBeskin fc Nokhrinal (12006D that, although e//To 1 on 
the (outer) fast surface, the perturbation method is not 
applicable beyond the fast point. As a result, we can 
only discuss the flow solution up to the outer fast point. 

The following assumptions for a PFD GRMHD flow 
along a hole-threading field line are considered. First, 

2 Due to a similar toroidal surface current distribu- 
tion, I = C/Aizr, on the equatorial plane, the parabolic 
force-free field around a black hole considered in BZ77: 
^ = y {r(l — cos 6) + 2 (1 + cos 6) [1 — In (1 + cos 0)]}, also fol¬ 
lows '3/ oc r( 1 — cosQ ) at large distance. This is one of the solutions 
of the source-free Maxwell equation in Schwarzschild spacetime: 

( (1 7% T,) 4', r ), T . + ( r 2 sing 'Fa),fl = while Equation g2) is one of 
the solutions of the source-free Maxwell equation in flat spacetime: 

(sue' 1 ':’').’' + = o. 
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inflow outflow 



Fig. 2.- The toroidal field, B as a function of magnetization parameter, oy, near the separation surface, r s for the inflow (left) and 
outflow (right) part of the solution. The superscript and ”+” denote the value computed at r -> r s _ and r -4 r 4 , respectively. The 
matching condition constrain the outflow by the inflow; that is, singling out the outflow solution that satisfy (Bj,) 4 = ( ) — (see §3.2 and 
§4.2). When 5 = 1 ((oy) - = (o*) + ft 13700), the matching condition is always satisfied. 


we assume ry = r s . The assumption of ry = r s and 
u p|r* = 0 ensures u p has a smooth transition from u p < 0 
(inflow) to ttp >0 (outflow). Second, to guarantee 
the flow is PFD. we re quire ay ^ 1 (see also §3.4.1 of 
iFendt fe Greinen (120011) for an estimation) and e -C 1. 
Furthermore, we assume e is constant along field lines. 
The higher the value of oy, the more magnetically dom¬ 
inated the flows are. 

Among all the parameter space we seek for the parame¬ 
ter set {Op, 77 , E, E}, which gives a similar time-averaged 
GRMHD simulation result in M06 for comparison. We 
therefore focus on a spinning black hole with its dimen¬ 
sionless spin a = 0.9 and a field line that threads the 
event horizon at mid-latitude, 6 = 60°. As mentioned 
in §2.2, the set {Bp, 77, E, L} can be equivalently deter¬ 
mined by { flp, oy, ry, r a}■ We adopt r* = r s (assumed), 
and Op = 1/2 Oh (similar to the result of M06) in both 
inflow and outflow region. Then we determine ay and 7 ~a 
( note that once ?'a is determined, the location of the fast 
surface is determined accordingly) by the constraint of: 
(i) (-Eem/Afl) ~ 10 2 near the separation surface^ simi¬ 
lar to the case in M06, and (ii) the matching condition. 

We note that Op ft 1/2 Oh is self-consistently ob¬ 
tained in a steady PFD GRMHD jlow solution for a 
monopole field geometry dBeskin fc Kuznetsova! 12000 ). 
while it may not be relevant for a parabolic field geom¬ 
etry. BZ77 examined the parabolic streamline in which 
Op decreases when shifting the angle from close to the 
pole to equatorial plan. In M06, the field geometry be¬ 
comes almost monopolar in the vicinity of the horizon, 
so that Bp ft 1 /2 Bh is observed along the field line (see 
also iBeskinl 120091 ). In the present paper, although the 
parabolic field is prescribed as a global field geometry, 
we nevertheless adopt the constant value of 1/2 Bh in 
our fiducial solution for convenience. 

4.2. Matching the Inflow and Outflow 

3 In top panel of Figure 7 of M06 is ~ 10 3 (see 

M06 for definitions), which equivalent to -Eem/Efl ~ 10 2 in terms 
of the definition in this paper. Note that our definition of Eem is 

47r smaller than V /)'' 11 used in M06. The factor of 4 tt is absorbed 
into the definition of in M06. 


Consistent solutions for PFD inflow and outflow along 
a field line are obtained iteratively until the matching 
condition is satisfied. For each ay, because Dp (= Bh/ 2) 
and ry (determined by where K' 0 = 0 is along the field 
line; see Appendix) are known, we can solve WE (equa¬ 
tion (13411 ) for the flow solution by requiring the physical 
solution to pass through the fast surface. For example, 
for the case e = 0.065, the profiles of (B^) - and (B ( j > ) + 
as a function of (ay) - and (oy) + , respectively, are shown 
in Figure [2] 

A consistent inflow/outflow solution exists when a suit¬ 
able set (e, (oy) - , (oy) + ) is applied. As mentioned in §3, 
the tendency is for {B/)~ ft const, to result in multi¬ 
ple choices of (o)~ such that (B^) - = {B/) + is satis¬ 
fied. This leads to certain amount of freedom for choos¬ 
ing the value for S. For simplicity, 5 = 1 is assumed, 
so (ay) - = ay| + = ay. We can read from Figure [2] 
that (B^) - = {B ( j > ) + is satisfied when (e,ay) = (0.065, ft 
13700). 

By the same method, for any specific value of e (or ay), 
there is a corresponding ay (or e) that satisfies the match¬ 
ing condition. The quantitative relation shows that, as 
ay increases, e also increases. This implies that, because 
there is more mass loading onto the field, the field pro¬ 
gressively bunches up toward to the rotational axis of the 
black hole. Finally, after ay is chosen by the matching 
condition, the parameter set {Bp, ay, ry, ta} of the in¬ 
flow/outflow part of the solution is uniquely determined. 
The relaxation of the assumption 8 = 1 is discussed at 
the end of §4.3.2. 

4.3. Self-Consistent Inflow/Outflow Solution 
4.3.1. Flow Properties 

We adopt the above parameter set (e, oy) = (0.065, ft 
13700) as the fiducial model parameters, because 
the resulting flow solution satisfies our requirement 
(Bem/Bfl)^* ~ 10 2 (§4.1). The conserved quantities, 
{Sip, ? 7 , E, L}, of our fiducial flow solution are shown in 
Table 1. The mass loading 77 changes sign according to 
u r and (equation (HOD ) in inflow and outflow re¬ 
gions. Because the sign has no specific meaning, the 
absolute value \rj\ is shown. By the assumption 5 = 1, 
|(?7)inflow| — |(? 7 )outflow| (equation (|39|) ). 
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TABLE 1 

Conserved Quantities of the Fiducial Flow Solutions 


(e = 0.065, 

a=0.9 

cr* ~ 13700) ^ with S = 1 

Inflow Outflow 


~ 7 x 10“ 5 

E(V)/n 

~ -112 ~ 114 

L{*)/n 

~ -720 ~ 724 

n F W 

~ 0.157 


Notes. 

t a consistent inflow/outflow solution is obtained when a suitable 
set (e,cr*) is applied, such that the matching condition is satisfied 
(see §3 and §4.2). 


TABLE 2 

Properties of PFD GRMHD Flow Along the Same 
Hole-Threading Field Line 




Inflow Solution 

Outflow Solution 


u r 

< 0 

> 0 


u e 

> 0 

< 0 


u* 

> 0 

> 0 


u 1 

> 0 

> 0 

E = 

£fl + Eem 

< 0 

> 0 


Lfl 

> 0 

> 0 


Eem 

< 0 

> 0 

L = 

Lfl+ Lem 

< 0 

> 0 


Lfl 

> 0 

> 0 


Lem 

< 0 

> 0 

e r = 

- P r _i_ pr 
~ C FL "r C EM 

> 0 

> 0 


t fr 

C FL 

< 0 

> 0 


pr 

C EM 

> 0 

> 0 


Notes. 

t For a stationary GRMHD inflow solution along a hole-threading 
field line, £p L < 0 is satisfied. In contrast, for an inflow along a 
non hole-threading field line during transient phase, £p L > 0 is 
possible fe.g. IKoide et al.ll2002t[K o rn i ssa.ro y||2005f ■ 

In the inflow region ( u r < 0), both E < 0 and L < 0 
indicates that the energy and angular momentum of the 
black hole is extracted outward (£ r > 0 and C r > 0). 
E //i of the outflow gives the maximum possible value of 
the terminal Lorentz factor (equation (EUll l. Although 
(Bem)” = (-Eem) + under the assumption 5=1, the 
absolution value of E = Eym + -Efl for the inflow is 
slightly smaller than the value for the outflow. This is 
because, in the inflow region, the fluid component £fl 
( or Lfl) has an opposite sign with the electromagnetic 
component £em (or Lem), partly canceling the electro- 
magnetically extracted energy (or angular momentum); 
whereas in the outflow region, the fluid and the electro¬ 
magnetic components have the same sign, both carrying 
the energy and angular momentum outward. The gen¬ 
eral properties of different physical components for PFD 
inflows and outflows are provided in Table [5] 

The extraction of black hole rotation energy by the 
GRMHD inflow is also indicated by the location of the in¬ 
flow Alfven surface. A remarkable feature in GRMHD is 
the existence of a negative energy region: once the Alfven 
surface of an inflow resides inside such a region, the blac k 
hole energy is extracted outward (Takaha shi et al.l ll9901. 
The inner boundary of the negative energy region is the 
inner light surface, and the outer boundary is defined by 



Fig. 3. — Characteristic points of a fiducial PFD GRMHD in¬ 
flow. Toward the black hole: Alfven surfaces (filled cyan triangles), 
light surfaces (empty green circles), and fast surfaces (filled blue 
squares). The event horizon and static limit (the outer boundary 
of the ergosphere) are shown by the thin solid and dashed lines, re¬ 
spectively . Field lines are represented by the thick solid line. The 
Alfven surfaces are located inside the the negative energy region 
(shaded region), implying that the black hole energy is extracted 
outward. 

gtt+gt^F = 0. Thus, the region must be inside the ergo¬ 
sphere, where gtt > 0. As the flow becomes increasingly 
PFD, the location of the Alfven surface moves toward 
the light surface, finally entering the negative energy re¬ 
gion (see the Appendix). For PFD GRMHD inflow, the 
fast surface is located very close to, and almost coincides 
with, the black hole event horizon. This is why the PFD 
inflow solutions are all similar, as mentioned in §3. In 
Figure [3] we plot the locations of the Alfven surface and 
fast surface of the flow, which share the same features 
mentioned previously. 

4.3.2. Radial Structure 

Let us now show the fiducial flow solution up to the 
fast surface in Figure [4] and compared the result (espe¬ 
cially Figures 7 and 8) in M06. The top panel of Figure 
|4] shows the opening angle of the prescribed field, which 
roughly follows a single power law 9 oc r~ 0 ' , which is 

in general more collimated compare to the result of M06. 
The locations of the characteristic surfaces are overlap 
onto the profile. The Alfven surfaces are located close to 
the light surfaces, and the inner fast surface is located 
close to the horizon. Note that in M06 the opening angle 
has different slope at a different radial range (see Fig¬ 
ure 10 of M06). Instead, our prescribed field line follows 
a single power law. Nevertheless, with similar require¬ 
ments at the separation surface ((-Eem/-Efl)|t-* ~ 10 2 ), 
the fast surface of the outflow is located at several hun¬ 
dred r g from the black hole, which is similar to the result 
of M06. 

The second panel of Figure 0] shows the profiles of the 
electromagnetic energy component Eem and the fluid 
energy component Eyl (both normalized by g), and 
the Lorentz factor T. Inside the ergosphere, g tt < 0, 
r is ill-defined. Therefore, only the profile segments 
outside the ergosphere are plotted. At large distances, 
gtt —> —1, r — > —ut = Eyl/r- In addition, for a 
PFD flow, E ss Eym Eyl when launching, so the 
maximum possible value of the terminal Lorentz factor, 
Too = E/g ss Eym/R near the separation surface. As 




















H.-Y. Pu et al. 


a result, the profile of T along the streamline is there¬ 
fore related to the conversion from Eem to Lfl- In the 
acceleration region (> 50?' g ), T roughly follows oc r 0,6 , 
which is similar to the result of M06, and t he analyt¬ 
ical result of r oc r 0 5 obtained in dBeskin fo Nokhrinal 
120061 1. It is expected that a further acceleration is ex¬ 
pected to be take place beyo nd the fast surface due to 
the m agnetic nozzle effect fe.g. lCamenzindl[l989tlLi et al.l 
I1992H . The conversion efficiency from Poynting to kinetic 
energy, which can be approximated by T/Too, is closely 
related to the location of the fast surface. For example, 
when the fast surface is located at infinity, T/Poo ~ 0. 
For the outflow solution, P/Too < 0.1 up to the fast 
surface, which is located at ~ 300r g . It is also inter¬ 
esting to note that the flow has already reached modest 
Lorentz factors (r ~ 5) at the fast surface, and most 
of the Poynting energy has not yet been converted to ki¬ 
netic energy. Note that, despite the final value of T at the 
fast surface is similar to the result in M06, the Poynting 
energy in M06 at fast surface has already experienced a 
significant decay (more than one order of magnitude) up 
to the fast surface. The reason why the fluid energy is 
not correspondingly increasing may be due to dissipative 
processes. In the inner region beneath the separation 
surface, —ut = Lfl/m < 1, as expected because the fluid 
is strongly bounded by the black hole’s gravity. In the 
outer region beyond the separation surface, —ut > 1, 
which implyies that the fluid is unbound and an outflow 
occurs. 

Similar to the energy conversion between the fluid and 
the electromagnetic components, the increase of the fluid 
component of the angular momentum Lfl is at the ex¬ 
pense of the electromagnetic component of the angular 
momentum Lem ■ The profiles of Lfl and Lem (nor¬ 
malized by g) are shown in the third panel of Figure [4j 
Again, the profile of the fluid component Lye/ g = is 
consistent with result of M06, but the decreases of the 
Poynting component in the simulation are much larger 
than our semi-analytical solution. 

The radial and polar components of the four-velocity 
of the flow, u r and u e , can be calculated from equations 
m and (l23l) . with u p determined by the WE. The other 
two components of the four-velocity, u 4 and v?, can be 
obtained by solving 

- g(u t + u^SIf) = E - fl F L , (44) 

subject to the normalisation u a u a = — 1. The velocity 
components u r and u 6 change signs across r = r s , while 
the velocity components u^ and u* remain positive in 
both the inflow and outflow regions. The angular velocity 
of the fluid, f 1 = u^/v/, which follows the black hole’s 
rotation, is however always positive along the magnetic 
field line. At the separation surface, u r = u e = 0, and 
hence S2 = flp- 

The radial and toroidal components of the orthonormal 
velocity at large distance are given by 

u r = yfg7r u r , (45) 

u+ = y/g&u+ (46) 

as shown in the fourth panel of Figure U The profile 
of u r is quite similar to the result in M06, but u^ has 
a relatively steeper profile compare to the simulation re¬ 
sult. We suppose this is related to the field configuration 


beyond the fast point, where we are not able discuss in 
current prescribed field configuration. 

The orthonormal components of the magnetic fields at 
large distance can be defined by 

B r = ^~ r B r , (47) 

B* = ^B*. (48) 

Note that B r is given initially when solving the WE, and 
L? 0 , which is not initially known, can be determined af¬ 
ter solving the WE. The bottom panel of Figure [4] shows 
the profile of the pitch angle, tan - 1 (|LT/L?^|). Because 
B r and B^ are both functions of g u (see the Appendix), 
they quickly decrease and change sign when entering the 
ergosphere ( g t t > 0). As a result, \B r \ and \B^\ are 
ill-defined close to the black hole, and we only plot the 
profile in the region where gu < 0. The reason why the 
pitch angle profile in M06 does not have this problem 
should be related to the definition of the field. The ex¬ 
plicit form of the magnetic field we adopt is provided 
in the Appendix. Nevertheless, at far region (e.g. the 
outflow region), spacetime becomes more flat and the 
differences of the definition are less important, our re¬ 
sult agrees with the result of M06. The locations where 
\B r \ = \B^\ are close to the light surface. At a large 
distance, \B r \ is well-described by \B r \ ss \B^\Rl/ y/g^, 
where Re = 1 /Qf, as also obtained in M06. 

At the end of this section, we discuss how would the 
flow solution would change if we adopt a <5, which also 
satisfies the matching condition, but does not equal to 
unity. Keep in mind that the outflow solution is well 
constrained by the matching condition, and the uncer¬ 
tainty of 8 is due to the degeneracy of the inflow solu¬ 
tions (§3). As a result, the outflow solution will remain 
the same if a different value of 8 is adopted. For the PFD 
GRMHD flow, because the location of the Alfven surface 
is always located near the inner light surface and the fast 
surface is always located close to the horizon, the flow dy¬ 
namics will therefore be similar. That is, u r ,u e ,u^ > ,u t , 
and therefore Lel/V = — «t and Lye/ g = will re¬ 
main almost unchanged. In addition, B r (prescribed) 
and B$ (constrained by the Znajek’s condition on hori¬ 
zon described in §3) will also remain similar. The elec¬ 
tromagnetic component, Lem and Lem, due to the de¬ 
pendence of OC I/ 77 , follow | (L/EM(inflow/(L E M)outflow | = 
| (LEM)inflow/(LEM)outflow | — | (^(outflow / (^?)inflow | — 1/(5. 

5. SUMMARY 

A semi-analytical scheme is presented to investigate 
the cold, PFD GRMHD flow solution along a Kerr black 
hole-threading field. The continuity of the outward 
Poynting energy flux across the separation surface is used 
as the matching condition to connect the inflow and out¬ 
flow parts of a PFD GRMHD fl ow solutio n. We con sider 
the parabolic field line of lBeskin fc Nokhrinal (l2006ll , and 
therefore the resulting flow passes through all the critical 
points at a finite distance. 

With similar black hole spin, angular velocity of the 
field, and magnetization at the separation surface, we are 
able to obtain a specific parameter set {flp, cr+, r*, ta} 
that gives inflow and outflow solutions in agreement with 
the time-averaged flow properties along a mid-level field 
line reported in the GRMHD simulation of M06. 
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Fig. 4.— Fiducial PFD GRMHD flow solution properties along a field line. Top panel: Jet opening angle of the prescribed field. Location 
of characteristic surfaces are also shown: separation point (plus sign), light surfaces (empty circles), Alfven surfaces (filled triangles), and 
fast surfaces (filled squares). The thin vertical line indicates the angular profile of the static limit (gtt = 0). Second panel: Electromagnetic 
energy component (upper solid line), L?em> and fluid energy component (lower solid line), Ep l> of the total energy E = const. = -Eem + ^fl? 
in unit of fluid rest mass energy. The profile of F is shown only when gtt < 0 (dashed line). Third panel: Electromagnetic (upper solid line), 
Lemj and fluid (lower solid line), Lfl components of total angular momentum, L = const. = Lem + Lfl- Fourth panel: The orthonormal 
velocities u r and u^. Bottom panel: The pitch angle of the orthonormal field, tan -1 (\B r /B^\) (solid line), which is well-described by 
tan -1 \Rl / (dashed line) at large distance, where Rl — 1 /Ftp- Because the orthonomal field is related to gtt and becomes ill-defined 

near the black hole, the pitch angle is only shown when gtt < 0. Along the field line, the location of the event horizon, the static limit, and 
the separation point are indicated by the vertical solid, dotted-dashed, and dashed lines, respectively. 
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In this current work, due to the limitation of the pre¬ 
scribed field configuration, we can only discuss the flow 
solution up to the outer fast surface. As a future work, a 
better consideration of the field configuration could help 
to explore of the flow acceleration beyond the fast sur¬ 
face, where the major jet acceleration takes place. 

Compared to the GRMHD and the general rel¬ 
ativistic f orce-free electro dynamics (GRFFE) (e.g. 

[McKinney fc Naravanl l2007f) numerical simulation ap¬ 
proaches, the semi-analytical approach provides a com¬ 
plementary understanding of the relativistic jets, in the 
sense that the numerical dissipative process is absent, 
and that the fluid component is included. The stationary 
solution obtained by the scheme can also be provided as 
a reference of the time-averaged GRMHD jet behaviour 
in numerical simulations. 

APPENDIX 

NOTES ON THE MAGNETIC FIELD 

Here we present the explicit form of the magnetic field. The covariant magnetic field defined in Equation (1131) 

Bp = , (Al) 

can be alternatively written as 

B » = i^ Fa ^ v , (A2) 


where ej, Ma/ g = \J~-g [vgaf3\, and [vfj,ap\, with \J~^g = Esinfl. Since = (1,0,0,0) and = 

(9tt, 0,0, we can quickly read from the above definitions that B t = 0 but B l ^ 0. 

The components of the magnetic field are therefore given by 


Br = V=9F<>*, 

(A3) 

nr 9tt + ^F<? t(f> t-, 

D = ,- red) , 

(A4) 

Be = , 

(A5) 

D 0 9tt + 9t<l> r 

= /- r r( k , 

V=9 

(A6) 

and 

B^ = V=gF re , 

(A7) 

B t = -3^ Fr6 . 

V^g 

(AS) 

With the relations 

p H= g ds9tt + nF 2 9t * Fet, 

Asm 0 

(A9) 

TpT(j> n rr9tt 9t<f> j-p 

r ~ 9 A • 2/1 » 

A sin 0 

(A10) 

rprO rr 06 jp 

F =99 F re , 

(All) 

one can check B r = g rr B r , B 8 = g 88 B s , and B$ = g^Bf Note that, despite is 

ill-defined near a Kerr black hole because gtt changes sign when entering the ergosphere. 
At large distance, the metric become Minkowski spacetime in spherical coordinates, 

finite at all region, B M is 

ds 2 = -dt 2 + dr 2 + r 2 d0 2 + r 2 sin 2 Odcf 2 , 

(A12) 

and y/—g = r 2 sin0. In this limit, the orthonormal field has the form 


B — i/g^B — . Fet , 

r z sm 6 

(A13) 
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(A14) 

(A15) 


In the following we outline the characteristic surfaces of cold GRMHD flow, including the light surfaces, the separation 
surface, and the Alfven and fast surfaces. 


Light Surfaces 

The surfaces defined by Kq = 0 are the light surfaces. There are two light surfaces in a black-lrole magnetosphere, 
the outer and the inner light surfaces. In the regions outside the light surfaces (where Kq < 0) the fluid streams 
radially so as to avoid the toroidal velocity exceeding the speed of light. The outer light surface is formed in the same 
manner as the light cylinder in a pulsar magnetosphere, but it does not necessarily have a cylindrical shape in a Kerr 
space-time. The inner light surface is formed due to strong gravity. Only when the black hole and the field line are 
not rotating does the inner light surface coincide with the black hole event horizon. 


Separation Surface 

In the cold limit the fluid acceleration along a field line, u' p (where prime denotes the derivative along the flow 
streamli ne), changes direction at a certain point. The location, r s , at which the change occur forms a separation 
surface dTakahashi et al.lll990l: iHirotani et all(1992 ). The fluid, starting with negligible velocity at r s , is accelerated 
inward inside the separation surface, creating an inflow. It is however accelerated outward outside the surface and 
develops an outflow. 

The separation surface is inside the region bounded between the two light surfaces, and is determined via searching 
for where K 0 = 0 along each flow streamline in the calculations. Figure [5] shows how r s on a specific field line (flow 
streamline) is determined in the demonstrative case with Kq{cl, Of) = A'o(0.9, Oh/2), (where Oh is the angular velocity 
of the black hole and r + is the radius of the outer event horizon). The location where Kq = 0 and K' 0 = 0 along 
the field line can be read from the contours of Kq, which are part of the light surfaces and the separation surface, 
respectively. Note that the locations of the light surfaces and the separation surfaces are independent of the flow 
parameters, such as the mass loading, as they are determined only by Ao(a, Op) and its derivative, Kq. 


Critical Surfaces for cold GRMHD flows 

Critical points appear when D vanishes in the expression of (lii-Up/ = N/D. In the cold limit, there are two critical 
points. The Alfven critical point corresponds to where u p is equal to the poloidal Alfven speed, i.e. 


u iw( r ; 



Airgm 


An 


(Bl) 


and the fast magnetosonic critical point corresponds to where u p equals the fast magnetosonic speed, i.e. 


(see iTakahashi et al.1 fl990 ). 
At the Alfven surface 

Setting Up = u\ w yields 


4m (G v k) — 4w 


Bi 


4-7r/mA sur 6 


M\ = K 0 


47t// V ° V 


(B2) 


(B3) 

(B4) 


Since n and ?y are positive, Kq > 0 at the Alfven surface. The Alfven surfaces are therefore constrained inside the 
region bounded by the light surfaces (where Kq = 0). In addition, K 0 —> 0 as 77 —> 0, implying that the Alfven surfaces 
approach the light surfaces when mass loading decreases. 

Because the flow must be super-Alfvenic outside the light surfaces (when shocks are absent), would the flows 
downstream, outside of the light surfaces eventually reach fast magneto-sonic speeds? The answer to the above 
question is different for inflows and outflows. For the inflow, the magneto-sonic speed is certainly reached, as causality 
requires that th e flow speed must surp ass all the possible characteristic speeds before the flow would enter the black hole 
event horizon (Takahas hi et al.l|19 90). For the outflow, whe ther or not the flow speed w ill reach the fast magneto-sonic 
speed depends on how fast the field decays along the flow (ITakahashi fe Shi bata 1993). 

If the fast surface exists, the physical flow solution for the WE can be uniquely determined after specifying three of 
the conserved quantities, and searching for the last one until the the flow can smoothly pass the fast surface, (see, e.g. 
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0 5 10 

r sine / (GM/c 2 ) 



Fig. 5.— Determining the location of the light surfaces and the separation surface for a field line. Top : Contour plot of Kq for the case 
a = 0.9 and Qp — ^h/2. The black area represents the region enclosed by the black hole event horizon. The contour of Kq = 0 is indicated 
by the solid line and the contours in region Kq > 0 are indicated by the dotted lines, with values at Kq = 0.1, 0.2, ..., 0.8, 0.9. The 
thick solid line is a representative large-scale, black hole-threading field line. The two green circles on it indicate the locations of the light 
surfaces, which corresponds to Kq = 0 and the red cross indicates the location of the the separation surface at which Kq = 0. Bottom : 
The value of Kq along the field line. The locations of the light surfaces and separation surface are indicated by the same symbols. 


Appendix C in lPu et all (120121 1 for the case of inflow as a demonstration). At the fast surface, where u 2 = u| M , 
have 


we 


n 

47r/i 


Kn 


B l 


Bp A sin 2 9 ) 


= V 


(B5) 


B y equati on (fTUl) . while all else being equal, a relatively smaller r) is expected to produce a stronger B$ (oc F r g) (see 
iPudritz et al.l ( 2006!) for a Newtonian version of such MHD feature). As a result, a smaller Kq is required to satisfy 
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equation (IB5[) when a smaller mass loading is applied. That is, the location of the fast surface moves farther away 
from the light surface as the mass loading decreases. For a GRMHD inflow, the location of the fast surface gets more 
and more closer to the event horizon. 
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